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In the DAE5ALUS (Decay- At-rest Experiment for 5c p studies At the Laboratory for Underground 
Science) project, high power cyclotron chains are proposed to efficiently provide proton beams 
with a kinetic energy of 800 MeV and an average power in the MW range. Space charge plays a 
pivotal role in both the injector and the ring cyclotrons. Large-scale particle simulations show that 
the injector cyclotron is a space charge dominated cyclotron and that a 5mA beam current can be 
extracted with tolerable beam losses on the septum. In contrast, in the ring cyclotron, no space 
charge induced beam loss is observed during acceleration and extraction. 

PACS numbers: 29.20.Hm;29.27.Bd;41.20.Cv 



I. INTRODUCTION 

The DAEDALUS collaboration is designing advanced cyclotrons that accelerate molecular hydrogen ions to produce 
decay-at-rest neutrino beams for a novel search for CP violation in the neutrino sector [1] . The DAEDALUS experiment 
calls for measurements at three values: L = 1.5 km, 8 km and 20 km. Three accelerator modules, called near-, mid- 
and far-site are used, providing 800 MeV protons with 0.8, 1.6 and 4.8 MW average beam power respectively. The 
experiment imposes a strong constraint on each accelerator in the complex, that is, to operate with a 20% duty cycle 
that allows the identification of the source of a neutrino event, and the measurement of the cosmogenic backgrounds 
[2]. Consequently, the current during beam-on should be 5 times the average beam current; i.e., in the mid-site 
accelerator module, 2 mA of average current requires 10 mA of peak current. The additional constraints on the 
accelerators of compactness and cost-effectiveness are not addressed in this paper. 
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This paper describes the beam dynamics study of the mid-site cyclotrons. The machine configuration discussed here, 
shown schematically in Fig. 1, is the mid-site cyclotron complex consisting of two cascaded cyclotrons. The injector 
cyclotron (DIC) is a four-sector compact machine, which accelerates a beam of rlj up to 60 MeV/amu. The beam is 
then extracted by an electrostatic deflector and is transported and injected into an eight-sector superconducting ring 
cyclotron (DSRC), in which the beam is accelerated to 800 MeV/amu by four single-gap rf-cavities. Two stripper 
foils can be used to extract two proton beams at the same time from the ring cyclotron, however, in Fig. 1 for sake of 
simplicity just one trajectory is shown. The key parameters of these two cyclotrons are listed in Table I. The general 
characteristics of all the sub-systems are elaborated in the Ref. [3-5] . 



# t Compact Injector 




FIG. 1: A schematic overview of the accelerator complex at the mid-site accelerator module, which is located at 8 km away 
from the underground detector. Only one extraction channel from the DSRC is shown in this schematic drawing. 

One of the main challenges of this high power accelerator complex concerns the control of space charge effects and 
the beam loss rates at extraction for both cyclotrons. The choice of accelerating instead of protons has two major 
advantages: a) The extraction of in the DSRC with a stripping foil does not require a clean turn separation at 
extraction, which would be mandatory for bare protons, b) In addition the space charge effects for are smaller. 
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TABLE I: Key parameters of the mid-site cyclotrons 





orbit radius 


kin. energy 


ctvg. power 


avg. field harmonic sector cavity 


cavity turn 




(cm) 


(MeV/amu) 


(MW) 


(T) 


no. type 


no. no. 


DIC normal 


5. . . 200 


0.035. ..60 


0.12 


0.95. ..1.17 6 


4 double-gap 


4 107 


DSRC superconducting 


190. ..480 


60. . . 800 


1.6 


1.06. ..1.88 6 


8 single-gap 


4 401 



A measure of the strength of space charge is the generalized perveance [6] 

K= qI 3R3 , (1) 
27T£om7 J 

where q, I, to, 7 and (3 are respectively the charge, current, rest mass and relativistic parameters of the particle beam. 
The higher the value of K, the stronger the space charge effects are. According to Eq.(l) the K factor for 5 mA, 800 
MeV/amu beam in the DSRC is equal to that of 2.5 mA, 800 MeV proton beam with the same 7. In consequence 
it is similar to the space charge effects that are present in the 2.4 mA, 590 MeV proton ring cyclotron at PSI. Similarly, 
the perveance of a 5 mA, 35 keV/amu H^" in the DIC is approximately equal to that of a 2 mA proton beam at 30 
keV, which has already been achieved in commercial compact cyclotrons, such as CYCLONE-30 [7] and TR-30 [8]. 
Therefore, a 5 mA beam seems not far beyond contemporary cyclotron capabilities with respect to space charge. 

This simple analysis, however, offers only the simplest understanding about space charge effects. Many other 
factors, such as the external focusing, the path length in cyclotrons, have significant impacts on the overall effects 
of space charge and the associated beam dynamics. Therefore, one can not get a complete picture without precise 
calculations of beam dynamics. A practical project demands detailed prediction of the effects of space charge on the 
beam's evolution based on accurate beam dynamics simulations of these high power cyclotrons. The stringent systems 
requirements of limiting controlled and uncontrolled particles losses requires simulations and measurements with more 
than three to four orders of magnitude in dynamic range. Simulation and modeling using large scale optimization 
techniques will be future research directions to meet the stringent requirements with respect to precise models and 
online control. 

The space charge effects in the DAEDALUS cyclotrons have been studied quantitatively by numerical methods 
implemented in the parallel code OPAL [9]. The beam dynamics model used in the simulation is described in 
Ref. [10, 11]. The report of the first design [5] has been followed by several iterations of magnet design and space 
charge simulations. Accordingly, the cyclotron structures have been improved to reduce beam loss to the acceptable 
range. The following sections discuss the space charge simulations based on the latest design solution of both DIC 
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and DSRC [3]. 

II. SPACE CHARGE IN THE DIC CYCLOTRON 

The maximum beam power delivered by the injector is determined mainly by space charge effects. Space charge 
forces can increase the beam emittance and consequently reduce the extraction efficiency. Reduced efficiency could be 
a serious limitation because beam losses must be kept below 200 W. This limit, following the experience of PSI and 
other laboratories, is found to be a practical maximum to allow for routine hands-on maintenance in the cyclotron 
vault. For the DIC, 200 W corresponds to an allowed beam loss of 1.67 /iA of H^~ beam at extraction. Hence one 
must pay close attention when designing for clean extraction from the injector cyclotron. 

A. Stationary matched beam formation 

As is described in literature [12-17], the space charge force combined with the radial-longitudinal coupling motion 
develops a vortex motion inside the bunch that can change the shape of the beam. If the vortex motion in the 
horizontal-longitudinal plane is strong enough, the beam will be approximately circular in this plane, as shown in 
[10, 18] and experimentally verified [19] at PSI. In the PSI Injector II cyclotron a stationary compact beam is developed 
within the first several turns, which remains essentially unchanged until extraction and the beam phase width is about 
2° at the last turn. Therefore, the two flattop cavities of this cyclotron, that were originally designed to suppress the 
increase of energy spread for the long beam, are now operating as accelerating cavities. With this setting, Injector II 
can deliver up to 3 mA of proton beam for the PSI Ring Cyclotron. To evaluate the necessity of the flattop cavity 
for the DIC cyclotron, one must investigate whether this working mode exists. 

When a coasting beam is perfectly matched with the cyclotron, its beam size oscillates periodically with a frequency 
equal to the number of sectors. Baumgarten [20] developed a theoretical model to compute the second moments of 
an ellipsoidal distribution (er matrix) matched to a cyclotron, in a linear approximation for a given beam current 
and emittance. In practice, the beam does not necessarily maintain an uniform distribution. Hence there exists non- 
linear components that we have to take into account. The full 3D code OPAL provides us with numerical methods to 
compute not only the a matrix, but also realistic particle distributions of a matched beam which includes nonlinear 
space charge forces. 

To search for a matched distribution, including non-linear space charge forces, we use as initial guess a spherical 



symmetrical particle distribution as suggested in Ref. [20] with the beam radius expressed by 

(To (l + f - fj) for < a < 5/2 



cr (1 + a) 



1/3 



(2) 



for 5/2 < a 



together with 



i / \ 3 / 2 

L -r{—} and oo = 



5V10m 



(3) 



where c is the speed of light, e„ the normalized emittance, a — c/ljq the cyclotron radius, / the beam current, h the 
harmonic number, and /j,o vacuum permeability. We start by tracking a coasting beam with 5 mA, 1.5 McV/amu 
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FIG. 2: The rms size snapshot at 0° azimuth in the 100 turns (left), and the rms envelop in the last two turns (right). For 
this simulation, 10 s particles and mesh of 32 3 grid points is used. 



until the equipartitioning process has resulted in a matched, stationary distribution. The tracking result shows that 
the sizes of the initial unmatched beam converge to almost constant values within approximately 50 turns. The 
longitudinal beam size approaches the horizontal size, as is shown in Fig. 2 (left). At different azimuths for the last 
two turns, the beam size oscillates with a frequency of four per turn (equal to the sector number of the cyclotron), as 
is shown in Fig. 2 (right). This behavior indicates that the stationary matched distribution has already developed. 
Figure 3 shows the horizontal and longitudinal profile of the stationary distribution. One concludes that a) this 
distribution is quite different from a uniform distribution, hence space charge has nondinear component; b) compared 
with a Gaussian distribution with the same rms parameters, this distribution has a sharper peak, thinner waist and 
longer tail. 

Since the space charge force is proportional to the charge intensity in the bunch, the size of the final matched beam 
is a function of the beam current. Fig. 4 compares the stationary beam size for the current from 1 mA to 20 mA that 
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FIG. 3: The longitudinal beam profile of the matched distribution and its Gaussion-fit with the same rms size 

shows that both the beam size and emittance of the formed stationary beam are gradually increase with increasing 
beam current. The beam size increases along with current both in the horizontal and longitudinal directions, however, 
the rate of growth is decreasing. The beam rms sizes can be well fitted by a quadratic polynomial 

cr = al 2 + bl + c (1 mA < I < 20 mA), (4) 

where the values of a, b, c are -0.0024, 0.12, 1.40 in the radial direction and -0.0026, 0.12, 1.55 in the longitudinal 
direction respectively. The beam emittance is more sensitive to the beam current than is the beam size. As the beam 
current increases from 5 mA to 10 mA, the horizontal and longitudinal emittances increase by approximately 50%, 
but the horizontal and longitudinal sizes increase only by approximately 20%. 




Current (mA) Current (mA) 

FIG. 4: The stationary beam rms size (left) and the normalized rms emittance (right) vs beam current 



From the simulation, one concludes that the stationary matched distribution with a short length can be formed for 
a 5 mA beam current. In consequence, no flattop cavity is required, as is the case for the PSI Injector II cyclotron. 
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Therefore, the four valleys are available for installing the accelerating cavities maximizing the energy gain per turn 
and hence to obtain a sufficiently large turn separation at extraction to achieve the required low beam losses. 

B. Space charge effects during acceleration 

The Versatile Ion Source (VIS) [21, 22] will be used to provide 35 keV/amu beam for the DIC cyclotron. The 
VIS is an off-resonance microwave discharge source operating at 2.45 GHz. It has demonstrated good CW proton 
currents (>35 mA) with a normalized emittance less than 0.2 7r mm-mrad and an extraction voltage which can be 
raised up to 70 kV. 

The spiral inflector and central region of the DIC is in the design phase, and an experimental test of the central 
region is in preparation. The simulations of the DIC described in this paper start at the exit of the central region. 
Considering that both the space charge effects in the injection line and the transverse-longitudinal coupling motion in 
the spiral inflector inevitably increase the emittance, the initial normalized emittance at the exit of central region is 
set to three times larger than that of ion source, i.e., 0.6 7r mm-mrad. The phase acceptance and initial energy spread 
are assumed to be 10° — 20° and 0.4 % respectively. In order to reduce the tail particles of the extracted beam, four 
collimators are placed at around 1.9 MeV/amu to cut off about 10 % of the halo particles. With the four collimators, 
the simulation is carried out for 1 mA, 5 mA and 10 mA current respectively. To achieve the high statistics for 
halo particles, 10 6 macro particles are used the in the simulations and 64 3 rectangular grid are used for space charge 
calculation. The Hj" beam needs 107 turns to reach the final energy. The phase slipping during the acceleration is 
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FIG. 5: RF phase slipping of the central particle in the beam 

plotted in Fig. 5, showing the phase slipping is within ±10° in the first 100 turns. During the final turns it increased 
by 45°, caused by the a sharp field drop, needed to increase the turn separation by passing the v r —\ resonance. The 



total turn separation consists of two parts: 

dR dR 



dR 



, (accel) + ——(precession). 
dn dn dn 



The contribution of acceleration can be expressed by 



dR. , AE 7 1 

— (accel) = R x — x — — - x — , 
dn hi 7 + 1 v£ 



(5) 



(6) 



where R denoting the radius, AE is the energy gain per turn, E is the total energy of the particles, and v r is the 
radial focusing frequency. The maximum turn separation of precession motion can be expressed by 

dR, 



dn 



(precession) = 2 x x c x sin7r(l — v r ), 



(7) 



where x c is a coherent betatron oscillation amplitude of beam. Figure 6 shows the turn separation obtained by particle 
tracking. The separation of the last turn is ~ 20 mm, in very good agreement with Eq. (5), i.e. the contribution of 
acceleration and precession are 13 mm and 7 mm respectively. 
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FIG. 6: Turn separation of the central particle in the beam. 

Figure 7 shows the rms beam size during acceleration. Again the beam sizes changes rapidly during the first five 
turns. This time the mismatch is mainly caused by the acceleration, which was not taken into account when finding 
the matched distribution of the coasting beam. After an particle repartitioning process driven by space charge force, 
the horizontal and longitudinal dimensions converge to nearly same size and only slowly increase at higher energies. 
A direct consequence is the reduction of the longitudinal phase width and, as a result, the suppression of the energy 
spread. During the last few turns the longitudinal beam size increases rapidly because the longitudinal space charge 
focusing is destroyed by the phase slipping. However, the radial beam size is unaffected by the phase slipping, a 
point that is important for extraction. The vertical beam size gradually decreases because of the increase of external 
focusing. The effect of the four collimators on the rms beam size is small as expected. 
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FIG. 7: The rms size snapshot at 0° azimuth the during acceleration with and without collimation for a 5 mA beam. 

C. Single-turn extraction from the DIC 

In OPAL, a radial profile probe element enables one to compute the radial beam profile at the location of the 
electrostatic deflector. Figure 8 (left) shows the radial beam profile at the position of the electrostatic deflector during 
the last four turns. One sees that the implemented collimation scheme reduces the charge intensity at the deflector 
position by a factor of four — a large improvement with respect to controlled beam losses at the deflector. Figure 
8 (right) shows the radial-vertical particle distribution at the deflector position during the last two turns. A very 
conservative septum thickness of 0.5 mm keeps the beam loss at extraction to less than 120 W. In the DAEJALUS 
experiment, the DIC will work at the a duty cycle of 20%, hence the beam power at the deflector is less than 24 W. 
The Isotope Decay-at-Rest (IsoDAR) [23] experiment is an intermediate step of DAEDALUS, in which the DIC will 
be used as a low energy high intensity driver. The IsoDAR electron antineutrino source requires the DIC to work at 
a duty cycle of 90% , in which case the beam power will be less than 108 W. Based on the experience with the PSI 
Injector II cyclotron, we conclude that the most crucial part of the DIC, the extraction is feasible. 

To evaluate the sensitivity of the beam losses to variation of injection conditions, we simulated scenarios with 
different phase acceptances in the central region and different beam currents. Figure 9 (left) shows the radial beam 
profile during the last four turns for 1 mA, 5 mA and 10 mA average beam current with initially 20° phase width. 
The radial beam sizes broaden with increased current thereby increasing the beam loss on the septum. An interesting 
phenomenon is that at the deflector septum position, for the 1 mA beam has more halo particles (i.e., the higher 
relative beam loss) than that of the 5 mA case. This behaviour is attributed to the fact, that the 1 mA beam is at 
the transition between the emittance-dominated and space-charge-dominated. Therefore beam distribution is more 
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FIG. 8: The radial profile of the last 4 turns at the center on the valley for 5mA current (left) and the the r-z projection of 
last 2 turns particle distribution for the case with collimation (right). The total macro particle number in a bunch is 10 6 . 



easily compromised by the phase slipping. Fig. 9 (right) shows the exponential increase of the power of beam loss on 
the septum versus beam current, which is well fitted by the exponential function 



P = 2 (exp(//2) - 1) . 



(8) 
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FIG. 9: Left: the radial profile of the last 4 turns at the center on the valley for the different beam current with the initial 
phase width of 20°; right: the resultant beam loss power on septum at 20% duty cycle and the exponential fitting curve. The 
total macro particle number in a bunch is 10 6 . 



Table II summarizes the beam loss for different beam currents and injection phases. At the level of 5 mA, the beam 
loss on the septum is not sensitive to the injection phase width thereby allowing the flexibility for the central region 
design. When the current increases to 10 mA, the beam loss on septum increases tenfold to a level which exceeds the 
tolerable level of loss. 
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TABLE II: The beam loss power on the septum for 20% duty cycle (uints: W). 
injection phase width 1 mA 5 mA 10 mA 
10° 6 24 240 

20° 14 22 316 

III. SPACE CHARGE IN THE DSRC CYCLOTRON 

In the DSRC, the beam loss at extraction is relaxed by using stripping extraction which has an efficiency close to 
100%, even if the beam orbits at the extraction radius are not well separated. Production of neutral hydrogen H is a 
significant problem in stripping H~ , but for the beam the fraction of Hq is expected to be substantially lower, as 
there is only one electron to remove, and this electron is much less tightly bound (2.7 eV) than an electron in neutral 
hydrogen (13 eV). The main challenge comes from two other aspects: a) the vertical focusing at higher energies, which 
is usually weak in a superconducting cyclotron; b) the relatively large energy spread of the resulting proton beam, 
due to the nature of multi-turn extraction. The following section investigates these two questions. 

A. Space charge effects during acceleration 

Since the layout of the transfer line between the DIC and the DSRC is not yet fixed, in the simulation we assume 
a typical gaussian distribution with conservative values for the normalized rms emittances of 0.9 7rmm-mrad in both 
the transverse directions, full phase width of ±10° and energy spread of 0.1%. Four single-gap cavities with a peak 
voltage of 1 MV accelerates the Hj" beam in 401 turns up to an energy 800 MeV/amu. Figure 10 shows the maximum 
vertical extent of the beam during acceleration for different current scenarios. Although the vertical extent of the 
beam is insensitive to space charge, but the beam extents gradually to 30 mm (full width) because of the weak external 
focusing, starting around turn 220. As the height of the magnet gap is 80 mm, no significant beam loss is expected. 
Figure 11 shows the longitudinal-radial configuration space of the final turn at extraction; the beam extends to ^20 
cm for different intensities. A stationary, compact shape will not develop, because the strength of the space charge 
force scales with I/7 2 . Hence the needed force for the vortex motion is too small. For the cases of 5 mA and 10 mA, 
the beam splits into three sub-bunches longitudinally that exemplify the longitudinal instability introduced in Ref. 
[17, 24]. 

For high intensity cyclotrons with small turn separation, single bunch space charge effects are not the only contri- 
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FIG. 10: The maximum beam extent in the vertical direction in the acceleration 
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FIG. 11: Top view of the final turn for 1 mA (left), 5 mA (middle) and 10 mA (right) beam current 

bution to beam dynamics. Along with the increase of beam current, the mutual interaction of neighboring bunches 
in the radial direction becomes ever more important. In the DSRC, the turn separation at injection is 25 mm and 
gradually decreases to 3 mm at extraction. Therefore, the interaction of neighboring bunches is of greatest importance 
at large radii where the bunches overlap and where the extraction foil is located. We have used the "Start-to-Stop" 
model [11] which is implemented in OPAL code to study these effects in detail. 

The simulation starts with a single bunch and automatically transfers to the multi-bunch tracking mode, when 
radially neighboring bunch effects become important. Figure 12 shows the projection of phase space onto the vertical 
(Z) - radial (R) plane for a mA (left) and a 5 mA (right) case. Beam halos extends vertically to ±20 mm which still 
clears the magnet sector and vacuum chamber. We conclude from the simulation that in the DSRC space charge and 
neighboring bunch effects introduce vertical beam halo, but the influence of those effects is at an acceptable value. 



13 




-0.02 -0.01 0.01 0.02 -0.02 -0.01 0.01 0.02 

Z [m] Z [m] 



FIG. 12: The projection of beam intensity onto the vertical & radial plane. 

B. stripping 

In the current design of the DSRC, the two strippers are located at 35° and 215° azimuth respectively (as are 
shown in Fig. 14). The stripper-element in OPAL samples the 6D-particle-distribution, and converts one into 
two protons. OPAL contains no physics processes for particle foil interaction and the creation of associated electrons, 
however dedicated programs are available for this set of problems. Figure 13 shows histograms of different beam 
quantities at the stripper foil. For this particular data set, the last 15 turns are contribute, exemplifying the need of 
a truly multi-bunch simulation. 

At extraction the full energy spread is about 1.2%, which is 12 times the energy spread of at injection. In the 
longitudinal direction, we see a long tail and the full phase width is 10°. The normalized rms emittances are e r =l.l 
7rmm-mrad and e z =1.4 7rmm-mrad, as compared with the injection emittance of 0.9 7rmm-mrad for both planes. 
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FIG. 13: Histograms of the proton distribution on the foil. 
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The particle distribution on the stripper foil is taken as the initial distribution for the calculation of the foil 
temperature and lifetime calculation. These results will be reported in a subsequent paper. 

C. Proton beam extraction 

Since the extracted protons have only half of the magnetic rigidity, the beam will be bend first towards the center 
of the machine and then towards the exit port. In this process the particles will experience complicated electrical- 
magnetic fields, crossing fringe fields and rf-cavities. Figure 14 shows the reference proton trajectory, together with 
two selected orbits for illustration purposes. The simulation shows that the vertical beam size extends into the 
hill gap in the region of the exit port. In order to compensate for this, an extra dipole field with 1.6 kG/cm is applied 
at the inner free space to strengthen the vertical focusing, as is shown in Fig. 14. 




-4 -2 2 4 

x (m) 

FIG. 14: The two proton extraction trajectories(red and black) and the strippers at 35° and 215° (magenta point), together 
with the first and last turn of the reference H J acceleration orbit (blue) and the shape of the extra dipole filed (green rectangle) . 
Other intermediate acceleration turns are not shown for clearance. 

The proton beam is tracked from the particle distribution on the stripper to the exit port and the extraction path. 
Figure 15 shows the beam envelope along the extraction trajectory. It can be seen that the beam is well-focused 
vertically with the help of the extra dipole field. The energy spread enlarges the longitudinal beam size, significantly, 
while the dispersion increases the horizontal beam envelop to 40 mm. We expect that the aperture of the extraction 
line towards the target will accept this beam without significant difficulties. 
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FIG. 15: The proton beam envelops along the extraction trajectory. 
IV. CONCLUSION AND DISCUSSION 

Given the tremendous impact of space charge on the beam's behavior in high intensity cyclotrons, precise beam 
dynamics simulation has been included at the very beginning of the conceptual design of the DAE<5ALUS cyclotrons. 
Based on the space charge study reported here, the current version of the cyclotron structure is well optimized with 
respect to the beam loss. 

The simulation shows that beam transport in the DIC cyclotron is space-charge-dominated allowing the formation 
the stationary matched distribution with short length. Therefore, no flattop cavity is required and the four valleys 
are available for installing the accelerating cavities. The H^" beam can be extracted with beam loss on septum less 
than 150 W for 100% duty cycle, which will not result in serious problems at extraction. 

Furthermore, in the DSRC the important beam properties at the stripper are dominated by the initial conditions. 
Space charge, and radially neighboring bunches, introduce vertical beam halo, but the influence is at acceptable levels. 
The multi-turn stripping scheme relaxes the constraint of beam loss at extraction. However, with the side effects of 
a larger emittance and larger energy spread of the extracted beam. Fortunately, the proton beam tracking along the 
extraction path shows the beam is well focused and therefore no significant beam loss is expected. 

Furthermore, other sources of beam losses which are not discussed in this paper. These phenomena include Lorentz 
dissociation of unquenched vibrational states of particles, collisions with residual gas and incomplete foil stripping. 
Studies of these aspects will follow in future papers. 
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